DUDLEY KNOX LIS. .A'.' - 
NAVAL POSTGRADUATr. 
MONTEREY, CALIF, 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




SINGLE-STATION ASSESSMENTS OF THE SYNOPTIC- 
SCALE FORCING ON THE 
MARINE ATMOSPHERIC BOUNDARY LAYER 

by 

John Patrick Gleason 
June 1982 

Thesis Advisor; K. L. Davidson 



Approved for public release; distribution unlimited. 

T204i,59 



SECUHITV CLAttiriCATlOM QW THIS fWt»< 



Dmtm Entmrmd) 



REPORT DOCUMENTATION PACE 


READ INSTRUCTfONS 
BEFORE COMPLETING FORM 


r nI^ont numbcn 


2. OOVT ACCetSlON NO. 


3. NCClNieHT*S C AT ALOG NUMaCn 


4. TiTLC 

Single-Station Assessments of the Synoptic-Scale 
Forcing on the Marine Atmospheric Boundary Layer 


». TV(»C OF REPOMT A PERIOD COVERED 

Master's Thesis 
June 1982 


t. FCNPONMING ONG. NCNONT NUMaCN 


7. AUTNONi'»> 




i. CONTNACT ON GNANT NOMSENf*) 


JOHN P. GLEASON 






t. ^CNCONNING ONOikMiZATlON NAMC ANO AOONCM 

Naval Postgraduate School 
Monterey, California 93940 


1 0 • AO G A A M £ L* Cm C N Ti AO JCC^. TASiC 

AACA A «0AK unit NUMBCAS 


1 1, CONT»*OLUNO O^^lCt MAMC AMO AOORKtS 

Naval Postgraduate School 
Monterey, California 93940 


June 1982 


IS. NUMRER OF PACES 

57 


\t MONlTo*IN(i AdiMfiv NAMfc A A60Me»i<«/ dUftwX 










IS*. OeCL ASSiPiCATlON/OOWNGNAOtNG 
SCNCOULE 


l«. OISTMlOUTlOM $TAT*M«MT f»l «#>!• 






Approved for public release; distribution unlimited. 


.7. oisT-lOUTlOH STATCMtHT fo/ •». M Ol.ek 30, II dllltml 


If. SU^^LCMCMTANY NOTCS 


1*. KEY WOROS (Ctmlinu* on f»wM aIO* II I4^illr Or Al»«k 


r) 


Single- Station assessment 
Single-station forecasting 
Vertical velocity 


Atmospheric Planetary Boundary Layer 
Adiabatic Method 
Kinematic Method 
Q- method 


20. ABSTNACT (Comitmf •• mldm H IMmnHiy mm¥^ 




Knowledge of the large-scale vertical velocity is required to predict the 
evolution of the atmospheric planetary boudary layer (APBL) . Since naval 
operations are often conducted in data sparse regions, single-station 
assessments of the vertical velocity are necessary. Data to evaluate such 
assessments were obtained from rawinsondes taken at San Nicolas Island, 
California. Vertical velocity estimates obtained by vertical integration of 

(CONTINUED) 



DO 1473 COITION OF I NOV •• I* 0B80LCTC 

’ ■’*** ^ S/H 0I0J-01*-«60I I ItCUUlTY CU AttIFiCATlOM OW TMI» P <kOt Dmf 

1 



^cu— tTV CL n»f 9m909m4- 



20. ABSTRACT (continuation). 

the moisture budget equation (Q-method) and by the adiabatic method were 
used to initialize an APBL 24-hour prediction model. RMS error statistics 
on predicted inversion height, potential temperature, and specific humidity 
were computed for the forecasts and compared to the RMS errors of a 
persistence forecast. Calculation of the vertical velocity by the Q-method 
showed the most promise. However, no single-station assessment method 
improved on the persistence forecasts. 
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ABSTRACT 



“ONreREy,„^““*;esc„oc. 



Knowledge of the large-scale vertical velocity is 
required to predict the evolution of the atmospheric 
planetary boundary layer (APBL) . Since naval operations are 
often conducted in data sparse regions, single-station 
assessments of the vertical velocity are necessary. Data to 
evaluate such assessments were obtained from rawinsondes 
taken at San Nicolas Island, California. Vertical velocity 
estiirates obtained by vertical integration of the moisture 
budget equation (Q-method) and by the adiabatic method were 
used to initialize an APBL 24-hour prediction model. RMS 

error statistics on predicted inversion height, potential 
temperature, and specific humidity were computed for the 
forecasts and compared to the RMS errors of a persistence 
forecast. Calculation of the vertical velocity by the 
Q-method showed the most promise. However, no 

single-station assessment method improved on the persistence 
forecasts. 
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INTROD OCTI ON 



Cne aspect of naval operations is battle group opera- 
tions in a fixed geographical area. This is currently 
referred to as MODLOC. Several scenarios lend themselves to 
MODLCC operations: the Viet-Nam era carrier based strike 

operations, pre-positioning for a projection of power, or 
simply a show of force. Training operations, as well, are 
often conducted in fixed areas to avoid hazards to commer- 
cial shipping and air traffic with weapons firings, or for 
privacy in operations. 

k second aspect is the need for self sufficiency of a 
deployed force. Any autonomous ability to make environ- 
mental prediction within the battle group reduces the 
requirements for shorebased support, decreases delivery 
time for perishable information and reduces the demand on 
communication assets. 

A third aspect is the effect of the atmosphere on the 
perfcrmance of weapons and sensors that depend on electro- 
magnetic and electro-optical wave propagation. Meteorologi- 
cal factors which affect such propagation can change 
significantly in four to six hours. The importance of APBL 
prediction to naval operations lies in: 
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• Prediction of the existence and evolution of electro- 
magnetic ducts. 

• Prediction of the inversion height, where it is known 
that optical propagation is degraded 'due to turbulence. 

This thesis is an examination of the accuracy of single- 
station assessments of the large-scale vertical velocity, 
and the effectiveness of these assessments for an existing 
APBL prediction model. Three methods of computing the ver- 
tical velocity will be discussed. Vertical velocities cal- 
culated by single-station methods will be used in an APBL 
model to ascertain the relative merits of these methods. 
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II. BACKGROUN D 



The single- station forecasting problem probably has been 
of interest as a military requirement since the origin of 



armed conflict. Oliver and Oliver (1945) put single-station 



assessments into a naval perspective: 



During the last few years, wartime conditions have made it 
necessary for isolated combat units to issue forecasts in 
regions where no network of meteorological stations could 
be available. Frequently the data from several stations 
or from reconnaissance planes are available, but in some 
regions the forecaster must rely only on surface and 
upper-air observations made at his own station. This is 
particularly true in the case of ships at sea. Hence, it 
IS important to develop proficiency at extracting informa- 
tion with limited aerological data.... He must further be 
able to distinguish from solutions which are internally 
consistent, and those which are inconsistent and therefore 
to be discarded. 



In APBL assessments the assumption of horizontal homo- 
geneity (one-dimensionality) is often used. All APBL pre- 
diction models require, as an externally known parameter, 
the large-scale vertical velocity which direcly affects 
changes in the height of the APBL. Furthermore, the height 
of the APBL affects other properties of the boundary layer 
such as the temperature and humidity. 

Eata for this examination were obtained from rawinsondes 
taken at San Nicolas Island (SNI), California. SNI data are 
sufficiently representative of a marine environment for this 
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Battalino et al. 



(1979) have stated that condi 



purpose. 

tions at SHI may be marine, continental, or mixed depending 
on the circulation patterns. A decision to use the SNI data 
was made, in part, because of the availability of similar 
data from 7andenberg (VHD) Air Force Base, California and 
Point Mugu (PMTC) , California, so that examination of the 
horizontal variability of the area and divergence derived 
verticaJ. velocity values could be accomplished. Fig. 1 is 
provided for geographic orientation. 

WW WW 




Figure 1. Southern California Coastal Waters. (After 
Battalino et al. 1979.) 
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A. THE PREDICTION MODEL 

An APBL prediction model in which the mean vertical 
motion is an external parameter is used to establish the 
affect of the derived values of vertical velocity on boun- 
dary layer evolution. The model used is based on the one 
presented by Stage and Businger (1981) and was coded for 
operational use by Brower (1982) . The predictive equations 
are: 



dq/dt = h”‘ (w* q* 


(0) ♦ We^q) 


(la) 


d0e/dt = h-» (w»0e* 


(0) + Re - Rb + We^9e) 


(1b) 


dh/dt = Wp ♦ We 




(1c) 



Where g is the total water content (liquid and vapor) , h is 
the inversion height, 9e is the equivalent potential temper- 
ature, w'q'(O) is the surface flux of moisture, w'9J (0) is 
the surface buoyancy flux, A is the difference between the 
variable value immediately above the APBL and its well mixed 
value in the APBL, Rj. is the rate of radiative heat gain per 
unit area at the cloud base, and is the rate of radiative 
heat loss per unit area near the cloud top. is the 

entrainment velocity and Wq is the large-scale vertical 
velocity at the inversion height. 
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The model atmosphere is two layered, composed of a tur- 
bulent lower layer in which equivalent potential tempera- 
ture, specific humidity and wind speed are constant with 
height, and a stable upper atmosphere with constant specific 
humidity and buoyancy gradients. The layers are separated 
at the inversion height by a zero-order discontinuity 
(jump). Fig. 2 depicts how a sample sounding is character- 
ized in the model atmosphere. 

Surface fluxes of moisture and bouyancy are diagnosed 
from bulk formulae (Lilly 1968). The liquid water profile 
is computed adiabatically and reduced by about 30 percent to 
agree with empirical results (Fairall et al. 1981.) Long 
wave radiative cooling is treated using the Stephan-Boltzman 
law, and short wave radiative warming is treated using the 
delta- Ed ding ton approximation (Fairall et al. 1981). 

The large-scale vertical velocity, Hq , is treated as an 
external parameter, and is constant over the model run 
(24-hr). The system is closed with the specification of the 
entrainment velocity which is accomplished by the assumption 
that dissipation of turbulent kinetic energy (TKE) is a 
fixed fraction of the production of TKE (Stage and Businger 
1981). A 30-minute time step is used. 



14 



SEP 1901 




Figure 2. Sample sounding in the model atmosphere. Solid 

line represents data, dashed line represents the 
model atmosphere. 
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. MET HODS 

At least three methods which can be used to compute the 
large-scale vertical velocity are: 

• The kinematic method. 

• The adiabatic method. 

• Integration of the moisture budget equation. 

The kinematic method requires wind sounding data from at 
least three appropriately spaced stations. The adiabatic 
method and the integration of the moisture budget equation 
are truly single-station assessments. 

A. THE KINEMATIC METHOD 

The kinematic method can be used to compute the vertical 
velocity when nearly simultaneous (within two hours) sound- 
ings are available at three or more appropriately spaced 
stations such as VND, PMTC, and SNI (Fig. 1). Following 
Saucier (1955), the vertical velocity can be derived: 

^w/dz = - (3u/^x + ^v/^y) , (2) 

where u, v and w are the zonal, meridional and vertical wind 
components. The horizontal wind components can be linearly 
approximated by a Taylor series expansion as: 



+ (^U/dx) 


X + Ou/^y) y. 


(3a) 


+ Ov/^X) 


X ♦ (av/^y) y. 


(3b) 
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Eg. (3) suggests that the horizontal derivatives may be 
obtained from a regression of the horizontal wind components 
onto the horizontal plane (Fairall et al. 1981) . Decomposi- 
tion of the wind vector from the sounding data at the three 
stations (SNI, 7ND, PMTC) into zonal and meridional compo- 
nents is required. Vertical profiles of the horizontal wind 
components were constructed at a 200 meter vertical incre- 
ment by linear interpolation from the irregularly spaced 
data levels to those levels that are an integer multiple of 
200 meters. Hultiple linear regression of each component of 
the horizontal wind onto the horizontal plane is performed 
as f Clio ws ; 



u = a + bx ♦ cy. 


(4a) 


V = d ♦ ex fy. 


(4b) 



where a through f are the regression coefficients. 

Comparison of Eg. (3) with Eg. (4) indicates that the 
coefficients b and f give the zonal derivative of zonal wind 
and the meridional derivative of meridional wind at the cen- 
ter of a triangle formed by the three locations and at each 
vertical level where the regression is done. 

Vertical integration of Eg. (2) is performed to give a 
vertical profile of vertical velocity to three kilometres; 

3 km 

W(2) = - ♦ (f) ) AZ r (5) 
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where w(z=0) = 0 and Az - 200m. 



B. THE ADIABATIC METHOD 



For isentropic motions there is no time change of poten- 
tial temperature following the motion. Based on this 
assumption Saucier (1955) specifies the adiabatic method for 
estimating vertical velocity as: 



(59/dt) » V (V9) 
(^9/^z) 



( 6 ) 



where 8 is the potential temperature. Saucier further 
states that the local temperature change due to advection 
may be on the order of that due to vertical motion. In this 
study, calculations of vertical velocity are made both 
including and neglecting the advective term. Description of 
the advection calculation is left to a subsequent section. 

The vertical gradient is determined by a linear fir to 
the lapse rate of potential temperature above the inversion. 
The local time change of temperature is obtained at fixed 
heights (z = 1.0, 1.25, 1.5 km) from sequential soundings at 
the same location. These heights are above the inversion. 
It is assumed that the vertical velocity at the surface is 
zero. The vertical velocity at the inversion is taken as 
the average of the vertical velocity at the fixed heights 
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linearly scaled to the inversion height. Thus, a single- 
station assessment of the vertical velocity is obtained. 

C. INTEGRATION OF THE MOISTURE BUDGET 

Integration of the moisture budget equation to compute 
the vertical velocity was performed by Lenschow (1973). 
This method will hereafter be referred to as the Q-method, 
and is based on the assumption of well-mixed specific humid- 
ity in the boundary layer, and that changes occur due to 
fluxes at the sea surface and inversion only. 

Vertical integration of the moisture budget equation 
through the boundary layer and across the interface between 
boundary layer and the stable atmosphere above yields the 
follcwing expression for vertical velocity, w, at the level 
of the inversion, h; 

A(q+1) (dh/dt) - h(d (q+l)/dt) ♦ w*q« (0) 

w = , (7) 

^(q^l) + (^h/2) 

where q and 1 are vapor and liquid water contents, w'q' is 
the moisture flux, /S is the vertical gradient of (q+1) in 
the mixed layer, and A{q*l) is the difference between total 
moisture in the mixed layer and total moisture immediately 
above the mixed layer. with no liquid water, z «h, and 
neglecting moisture advection Eq. (7) is simplified to: 
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w 



^q(dh/dt) - h(dq/dt) ♦ w'g»(0) 



( 8 ) 



Lilly (1968) assumes that is zero. Lenschow (1973) found 
ySh/2 to be small and certainly much less than the moisture 
jump. 

Eq. (8) provides another independent assessment of the 
vertical velocity at the inversion. The time derivatives 
(dh/dt and dq/dt) are determined from sequential soundings. 
The moisture jump and inversion height are taken as averages 
for the two soundings, and the surface moisture flux is det- 
ermined by the bulk method, i.e. , 

= ChV^o(3s - qio) (9) 

The exchange coefficient, c^, is computed from the wind 
speed dependent relation of Large and Pond (1980). The sur- 
face specific humidity, q^, is determined via Tetens* for- 
mula for the vapor pressure (Buck 1981). 

D. DETERHIN&TION OF ADVECTION FROM THE THERMAL HIND 

Under the geostrophic assumption, the horizontal temper- 
ature gradient is directly proportional to the vertical 
shear of geostrophic wind. Following Byers (1974): 

V = f-Md<^/dx)p ; u = -f-‘ (d<^/^y)p » (10) 
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where v and u are the meridional and zonal components of the 
geostrophic wind, f is the Coriolis parameter, <i> is the 
geopctential and the horizontal derivatives are on the con- 
stant pressure surface. Vertical differentiation of 
Eg. (10) and substitution from the hydrostatic equation 
yields: 

^v/Sp = f-‘ (-a ) , (11a) 

^u/Sp = f-» ( Of ) , (11b) 

where a is the specific volume. Substitution from the 
equation of state and rearranging terms, the components of 
the temperature gradient are: 

(12a) 

(aT/^y)p = (fp/R ) (au/dp) , (12b) 

where R is the gas constant for dry "air. 

With a wind sounding, assumed to be geostrophic above 
the boundary layer, the horizontal temperature gradient can 
be specified and a single- static n assessment of the thermal 
advection term in Eq. (6) is obtained: - (V • VT) =- (V • V9) , at 

a constant level. 

This calculation is performed by subtracting the wind at 
the first level above the inversion from the wind at an 
upper level. The advecting wind is the mean wind in the 
layer. 



(aT/^x)p = -(fp/R ) (^v/ap) , 
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E, SINGLE-STATION ASSESSMENT SCORES 

The separate vertical velocities computed from single- 
station assessments by the adiabatic method and the Q-method 
were both used as the large-scale vertical velocity in the 
AP BL mod el . 

Eoth methods required use of finite differences in time 
to approximate time derivatives. The time increment used 
was the time difference between subsequent soundings at SNI; 
it varied from four to eighteen hours. To establish the 
vertical velocity for a model run, the time weighted average 
vertical velocity was used, as computed from the SNI data 
for the 24 hour period following model initialization. 

The model requires initial values of inversion height, 
and initial profiles of potential temperature and specific 
humidity. Sea surface temperature, required for computing 
surface moisture flux, was constant through the model day. 
Wind speed, required for computing surface moisture and 
bouyancy fluxes, was input as a linearly changing function 
of time. 

Three model runs were made on each set of initial condi- 
tions; one with vertical velocity by the Q-method, a second 
with vertical velocity by the adiabatic method and the last 
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with vertical velocity set to zero. The advective term of 
Eq. (6) it was found to be much less than the local time 
change of potential temperature. Thus, it was neglected. 

Boot Mean Square (RMS) error was computed for predicted 
inversion height, h, potential temperature, 0, and specific 
humidity, q, 

! n 

i=i 

where N is the number of observations, Xj is the observed 
value of a variable, and y, is the predicted value of a 
variable. Observed values were obtained from rawinsonde 
data available within the 24 hours following the time of the 
initial data. Typically, two or three soundings were avail- 
able for each 24- hour period following initialization. 
Predicted values were taken from the model output for the 
time CO rr.es pen ding to that of the observed values (± 15 

minutes, due to the 30 minute model time step). 
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IV . s YN OPT IC SITO ATION 



Two periods (cases) were considered. Case I was from 31 
August to 4 September 1981 and consisted of the analyses of 
nine model predictions verified against twenty observations. 
Case II was from 15 to 18 September 1981 and consisted of 
the analyses of eight predictions verified against seventeen 
observations. 

A. CASE I (31 AOG - 4 SEPT 1981) 

Throughout the period the eastern North Pacific Ocean 
was under a high pressure system that was moving slowly 
eastward. Southern California was dominated by a thermal 
trough^ leading to weak offshore flow at the surface, while 
California was to the south of the strongest westerlies at 
upper levels. Pigs. 3 through 6 show the National Meteoro- 
logical Center (NMC) surface and 500 millibar (mb) analyses 
for Case I. 

Selected prints of the GOES WEST Infra-red (IE) imagery 
indicate low broken to overcast cloudiness in the vicinity 
of SNI (Fig. 7). Surface observations at SNI indicated 
overcast conditions with fog at night, and scattered low 
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cloudiness with haze persisting after the late morning 
burn-off of the fog. 

A Pacific high pressure system dominated the eastern 
ocean area west of California on 31 August (Fig. 3). A low 
pressure system was located in the Gulf of Alaska with its 
associated weak frontal system. Tropical Storm Irwin was 
dissipating southwest of the Baja Peninsula. A thermal 
trough dominated the central California valley. Very weak 
westerly flow aloft was predominant over southern California 
and most of the eastern Pacific Ocean. 

After 31 August the Pacific high intensified^ while the 
low pressure system in the Gulf of Alaska moved into the 
northwest United States and the assosiated cold front dissi- 
pated. The thermal trough over inland California moved 
slightly eastward and dominated the high desert region of 
Arizona and southeastern California (Pigs. 4 through 6). 
Flow aloft remained very weak over southern California and 
the low level flow in the vicinity of SNI was very weak (3 
to 5 kt) and predominantly offshore. 
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Figure. 3. Surface and 500 mb analyses, 1200 GMT 31 Aug 1981. 
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Same as 
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Figure 5. Same as 3 except 1200 GMT 3 Sept 1981. 



>4 
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Figure 6. Same as 3 except 1200 GHT 4 Sept 1981 




1245 GMT 31 AUG 31 



1315 GMT 1 SEP 31 




1245 GMT 2 SEP 31 1315 GMT 3 SEP 31 

Figure 7. GOES BEST IB, CASE I. 
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B. CASE II (15-18 SEPT 1981) 

Throughout the period the polar front was located in the 
eastern part of the North Pacific Ocean and Gulf of Alaska. 
A high pressure system dominated the area southeast of the 
polar front. A thermal trough was located over the inland 
California valley and through much of the southwestern 
United States and Mexico. Figs. 9 through 13 show the NMC 
surface and 500 mb analyses for Case II. 

Satellite imagery indicates considerable clear area in 
the vicinity of SNI. Surface observations at SNI reported 
scattered low and high clouds with occasional reduced visi- 
bility due to fog and haze. Selected prints of the GOES 
WEST IR imagery are shown in Fig . 14 . 

Cn 1 4 September high pressure prevailed over the eastern 
Pacific Ocean. The thermal trough dominated the southwest- 
ern United States and the Mexican Plateau. A low pressure 
system with a weak frontal system existed in the Gulf of 
Alaska. Weak unorganized flow existed south of the strong- 
est westerlies in a large region of very flat height gradi- 
ents over most of California and the eastern Pacific Ocean 
(Fig- 9) . 
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Ey 1 5 September a large amplitude ridge aloft had built 

r 

over the western United States and Canada (Fig. 10) . The 
thermal trough weakened in central California as the ridge 
axis in the eastern Pacific Ocean migrated eastward. Low 
level flow in southern California was offshore and stronger 
(12 to 2 2 kt) than in Case I. 

C. DIFFERENCES BETWEEN CASES I AND II 



Comparison of the mixed layer and overlying data clearly 
indicates that different synoptic situations existed for 
the two cases. These differences are summarized in Table I 



and illustrated in Fig. 8 : 



TABLE I 

Comparison of the two cases 



mean h 
range of h 
mean 9 
mean A9 
mean q 
mean Aq 

mean 30/dz above h 
mean ^q/^z above h 



Case I 
604 m 

511 to 66 7 m 

16.2 C 

15.0 C 
8.3 g/k g 

-6.0 g/kg 

4.2 C/km 
-.35 g/(kg-km) 



Case II 
262 m 

145 to 357 m 

18.1 C 

14.1 C 

9.8 g/kg 
-4.4 g/kg 

2.8 C/km 
-.83 g/(kg-km) 



On the average, inversion height was nearly 2.5 times 
greater in Case I than in Case II. Case II had a warmer and 
more moist boundary layer with smaller magnitudes of jumps 
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in bcth potential temperature and specific humidity. The 
vertical gradient of potential temperature (stability) above 
the inversion is 1.5 times greater in Case I than Case II, 
and the vertical moisture gradient is nearly 2-5 times 
larger in case I than Case II. 



E 
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Figure 8. Average Model Atmosphere: Case I vs. Case II 

^ Solid: Case I, dash: Case II 
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and 500 mb analyses, 0000 GMT 14 Sept 1981. 
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Figure 11. Same as 9 except 0000 GMT 16 Sept 1981. 
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Figure 12. Same as 9 except 0000 GMT 17 Sept 1981. 
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Figure 13. Same as 9 except 0000 GMT 18 Sept 1981. 



1315 GMT 14 SEP 31 



1015 GMT 15 SEP 31 




1015 GMT 16 SEP 31 



0315 GMT 13 SEP 31 



Figure 14. GOES WEST IR, Case II. 
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RESULTS 



Direct measurement procedures do not exist to obtain 
vertical velocity values. Therefore an indirect evaluation 
procedure was nesessary. This was accomplished by using the 
RHS error (section III,E.) for model predictions of inver- 
sion height, potential temperature and specific humidity, 
with the large-scale vertical velocity specified by the 
adiatatic method, Q- method, vertical velocity equal to zero, 
and a persistence forecast. It is noted that the Q-method 
is based on measured values of the prediction variables 
(Eg. (1) ) . Therefore the model is expected to perform well 
with vertical velocities computed in this way. 

In both cases the inversion persisted for several days. 
Eq. (1c) illustrates that H acts to keep a lid on the boun- 
dary layer only when subsidence (W <0) is occurring since H 
provides only for increasing the inversion height, i. e., 
"detrainment” is not allowed. Then from a pragmatic view, 
negative vertical velocities (subsidence) are sought while 
positive ones are not. Large-scale subsidence is also 
expected, from a synoptic viewpoint in both cases, since SNI 
is under the summertime subtropical ridge. 
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Figure 15. Frequency of occurrence of subsidence. Case Ij solid: 

Q-method; dash: adiabatic method; double dash: kinematic 

method . 
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Table II gives the averages of the vertical velocities 
for each case calculated by the Q-method, the adiabatic 
method, and the kinematic method. The standard deviation is 
given in parentheses. Figs. 15 and 16 show the frequency of 
occurrence of calculated vertical velocity for Case I and 
Case II, respectively. assuming that subsidence existed at 
the inversion height, the adiabatic and kinematic methods 
are unable to diagnose subsidence. The Q-method consis- 
tently gives subsidence where subsidence is expected. Dis- 
cussion of possible sources of errors in these methods is 
left to a subsequent section. 



Case I 
Case H 



TABLE II 



Vertical velocity (cm/sec) comparisons 
Q-method Adiabatic method Kinematic method 




A. DISCUSSION OF CASE I 

The RMS errors for Case I appear in Table III . The 
average time from initialization to verification for this 
case was 15.6 hours. Persistence forecasts yield the least 
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BHS error for all variables (h, 9, and q) . 



When the model 



is initialized with the large-scale vertical velocity com- 
puted by the adiabatic method or when the vertical velocity 
is assumed to be zero, very poor results are obtained in the 
inversion height prediction. Vertical velocities computed 
by the Q-method result in better prediction of inversion 
height than the adiabatic method. 

TABLE III 
RHS Error ( Case I) 

h 9 q 

Persistence 61m 2.2C .6 g/kg 

W=0 167 m 2.2 C .8 g/kg 

Q-method 106 m 2.4 C ,9 g/kg 

Adiabatic method 249 m 2.2 C .7 g/kg 

The large errors in the inversion height for the W=0 
prediction occur because the boundary layer is constrained 
to growth by entrainment. The poor results of the adiabatic 
method are also to be expected because of two basic problems 
with this method. First, examination of Eq. (6) reveals 
that the lapse rate is assumed to be constant and the numer- 
ator is approximated by a finite difference in time over six 
to twelve hours. An uncertainty of ± one degree in poten- 
tial temperature can lead to an uncertainty in vertical vel- 
ocity of ±.5 (cm/sec) , assuming nominal values of the free 

* 

atmosphere lapse rate (2.5 to 5 C/Km), Since the actual 
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vertical velocity may be positive or negative and often on 
the crder .5 (cm/sec), the adiabatic method is not adequate 
to specify either the magnitude or the sign of the vertical 
velocity. Second, one must choose an appropriate time 
interval over which the finite difference applies. In prac- 
tice the time interval was that of the rawinsonde schedule; 
however, the soundings must be scheduled suffiently close in 
time to reveal changes in sign of the tendencies of the 
boundary layer variables. 

The Q-method also suffers from the disadvantages of 
finite difference in the specific humidity (dq/dt) and the 
inversion height (dh/dt) . However, the moisture flux term 
has no such problem. In addition, specification of the 
moisture and inversion height variation with time required 
for these calculations does not require profiling devices 
such as radiosondes, but may be obtained from surface 
instrumentation. 

B. CISCOSSION OF CASE II 

Average time from initialization to verification was 
14.6 hours. Table IV summarizes the results. 

None of the predictions did particularly well with 
respect to RMS error of inversion height. The same ordering 
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TABLE IV 

RMS Error (Case II) 



Adiabatic method 



Persistence 



11=0 

Q-method 



h' 

105 m 
375 m 
142 m 





of RMS errors for inversion height prediction occurs with 
the vertical velocity specification methods in this case. A 
major difference in these results is that the model pre- 
diction of specific humidity with vertical velocity by the 
Q-method exhibits a 24t error, much greater than any other. 
This is to be expected since the average inversion height 
was 262m, over 40* lower than in Case I, and in Eg. (la) the 
specific humidity tendency is inversely proportional to 
inversion height. Hence, low initial inversion heights, as 
was the case, coupled with predictions of inversion heights, 
which are too low, contribute to the over prediction of spe- 
cific humidity- The Q-method yields the least RMS error in 
potential temperature predictions. 

C. CISCDSSICN OF THE KINEMATIC METHOD 

Vertical velocities calculated by the kinematic method 
were positive more often than not as indicated in Pigs. 15 
and 16 . They were much larger in magnitude than those cal- 
culated by either the adiabatic method or the Q-method. 
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Therefore application of these to the model would have pro- 
duced larger RMS errors than any presented. Two reasons for 
the poor results are: 

• Improper specification of the bottom boundary condi- 
tion. 

• Small sample size for the regression approach. 

ever the ocean the bottom boundary condition for Eg. (2) 
is that the vertical velocity equals zero at the sea sur- 
face. This is a good assumption since the surface is flat 
with respect to a local earth tangent plane. However# since 
SNI, PMTC, and VND are coastal stations, the boundary condi- 
tion for Eg. (2) requires specification of the terrain 
(Dutton 1976). 

D. CONCLOSION 

Since both Case I and Case II represent small samples of 
observations the error statistics have little stability. 
However, they do suggest that vertical velocity calculated 
by the Q-method has the most merit as a single-station 
assessment, clearly, the assumption of no vertical motion 
or vertical motion derived from the adiabatic method has 
questionable value in APBL assessments. Fig. 17 shows a 
model prediction with vertical velocity computed by the 
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Q-method. Fig. 18 is the same model run but with upward 
vertical velocity computed by the adiabatic method. Compar- 
ison of these figures shows the good agreement in q pre- 
diction indicated in Table III . On the other hand, the 
superior result of the Q-method is well illustrated by com- 
parison of the evolution of the inversion height. The 
importance of diagnosing subsidence when it exists is evi- 
dent by this comparison. 

The assumption that SNI represents a marine environment 
is a possible source of error in the Q-method. Also, sea 
surface temperature has an indirect effect on the calcula- 
tion of vertical velocity through Eq. (9). The sea surface 
temperatures used were obtained from the Sea Surface Thermal 
Analyses of the National Oceanographic and Atmospheric 
Administration, rather than being in-situ measurements. 

Additional research is needed in specification of the 
large- scale forcing of the APBL. Vertical velocity by the 

Q-method appears promising, however larger samples must be 
obtained to gain statistical stability. This need not be an 
expensive venture. One only needs continuous measurements 
of sea surface temperature, air temperature, humidity and 
wind. That can be done on a near-continuous basis with 
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currently available instrumentation that does not place 



great demands on personnel. 
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Figure 18- Same as 17 except Adiabatic Method 



51 



LIST OF BEFERENCES 



Battalino,T. , J. Rosenthal, H. Hendon, 7.R. Noonkester, 
1979: ” Marine/Continental History of Aerosols at San 

Nicolas Island during CEHCOH-78 and OSP II", TP~79~33 . 
Pacific Hissile Test Center, Point Mugu, CA. 



Brower, D. A. , 1982; "Marine Atmospheric Boundary Layer and 
Inversion Forecast Model", NPS Master's Thesis, Naval 
Postgraduate School, Monterey, CA. , 13U pp. 



Buck, A. L. , 1981; "New Equations for Computing Vapor 
Pressure and Enhancement Factor". J. AppI. Meteor.. 20. 
pp. 1527-1532. 



Byers, H.H., 1974; 
McGraw-Hill, Inc., 



General Meteoroloov f4th Edition). 
New York, Ip: 175^^1 . 



Duttcn, J. A. , 
New York, pp. 









1976; 



McGraw-Hill, Inc., 



Fairall, C.H., K.L. Davidson, and G. E. Schacher, 1981: "A 
Review and Evaluation of Integrated Atmospheric Boundary- 
Layer Models", NPS-6 3-81-0 04 , Naval Postgraduate School, 
Monterey, CA. , 89 pp. 



Large, W.G. and S. Pond, 1980: "Open Ocean Mommentum Flux 

Measurements in Moderate to Strong Winds", J. Phys. 
Oceanoqr . , 11, pp. 324-331 . — — 



Lenschow, D. H., 1973; "Two Examples of Planetary Boundary 

Layer Modification over the Great Lakes", J. Atmos. Sci .. 
30, pp. 568-581. 



Lilly, D.K.. 1968: "Models of Cloud Topped Mixed Layers 

under a Strong Inversion", Quart. J, Roy. Meteor. Soc. , 94, 
pp. 292-309. ^ 



Oliver, V. J. and M,. B. 
Single Station Data", 
E. Bella y, N.R. Beers, 



Oliver, 1945: "Weather Analysis from 

Handbook of Meteorology . F. A. Berry, 
McGraw-Hiil, New York, pp. 858-879. 



Saucier, W.J., 1955: Principles of Meteorological Analysis . 

University of Chicago Press. 



Stage, S.A, and Joost A. Businger, 1981: "A Model for 

Entrainment into a Cloud Topped Marine Boundary Layer. Part 
I: Model Description and Application to a Cold-air Outbreak 

Episode", J. Atmos. Sci .. 38, pp . 2213-2229. 



52 



INITIAL DISTRIBOTION LIST 



No. Copies 



1. 



Defense Technical Information Center 
Cameron Station 
Alexandria, VA 22314 



2 



2. Superintendent 
Attn; 

Library, Code 0142 
Professor C.N.K Mooers, Code 68 
Professor R-J. Renard, Code 63 
Professor K.L. Davidson, code 63Ds 
Professor G.E. Schacher , Code 6lSq 
Professor J. Dyer, Code 61Dy 
Asst. Prof. R.w. Garwood, Code 68Gd 
Naval Postgraduate School 
Monterey, CA 93940 



2 

1 

1 

10 

10 

1 

1 



3. Director 1 

Naval Oceanography Division 
Naval Observatory 
34th and Massachusetts Ave. NW 
Washington, D.C. 20390 



4. Commander 1 

Naval Oceanography Command 
NSTL Station 
Bay St. Louis, MS 395 22 



5. Commanding Officer 1 

Naval Oceanographic Office 
NSTL Station 
Bay St. Louis, MS 39522 



6. Commanding Officer 1 

Fleet Numerical Oceanography Center 
Monterey, CA 93940 



7. Commanding Officer 1 

Naval Ocean Research and Development Activity 
NSTL station 
Bay St. Louis, MS 39522 



8. Commanding Officer 
Attn; 

Dr. A. Weinstein 1 

Dr. A. Goroch 1 

Dr. Steven Burke 1 

Mr. Sam Brand 1 

LT Mark Schultz, DSN 1 



Naval Environmental Prediction Research Facility 
Monterey, CA 93940 



53 



9 



1 



Chairman, Oceanography Department 
O.S. Naval Academy 
Anapolis, MD 21402 



10. Chief of Naval Besearch 1 

800 N. Quincy Street 

Arlington, VA 22217 

11. Office of Naval Research (Code 480) 1 

Naval Oceanography Research and Development Activity 
NSTL Station 

Bay St. Louis, MS 39522 

12. Scientific Liaison Office 1 

Office of Naval Research 

Scripps Institute of Oceanography 
La Jolla, CA 92037 

13. Library 1 

Scripps Institute of Oceanography 

P.O. 6ox 2367 
La Jolla, CA 92037 

14. Library 1 

Department of Oceanography 

University of Washington 
Seattle, WA 98105 

15. Library 1 

CICESE 

P.O. Box 4803 

San Ysidro, CA 92073 

16. Library 1 

School of Oceanograrphy 

Oregon State University 
Corvallis, OR 97331 

17. Commander 1 

Oceanographic Systems Pacific 

Box 1390 

Pearl Harbor, HI 96860 

18. Dr. C.W. Fairall 1 

BDM Corporation 

134 0 Munras St. 

Monterey, CA. 93940 

19. Mr. Don Spiel 1 

BDM Corporation 

134 0 Munras St. 

Monterey, CA 93940 

20. Dr. Paul Twitchell 1 

code 37 OC 

Naval Air Systems Command 
Washington, D.C. 20360 



54 



21 



1 



. Dr. Alex Shantla 
Code 3173 

Naval Weapons Center 
China Lake, CA 93555 



22- Dr- Barry Katz 1 

Cod e F4 2 

Naval Surface Weapons Center 
White Oak Laboratory 
Silver Springs, MD 20907 



23. Dr- J-H- Richter 1 

Code 53 2 

Naval Ocean Systems Center 
San Diego, CA 92152 

24- Dr. Lothar Ruhnke 1 

Code 8320 

Naval Research Laboratory 
Washington, D-C. 20375 

25. Mr. Herb Hitney 1 

Naval Ocean Systems Center 

Code 532 

San Diego, CA 92152 

26. Hr. Herb Hughes 1 

Naval Ocean Systems Center 

Code 532 

San Diego, CA 92152 

27- Hr- Stuart Gatham 1 

Code 83 20 

Naval Research Laboratory 
Washington, D-C. 20360 

28. Commander 1 

Naval Sea Systems Command 

PMS-4 05 

Washington, D.C. 20360 

29. Mr- Ted Zuba 1 

Code Air-370 

Naval Air Systems Command 
Washington, D-C. 20360 

30. Mr- J. Rosenthal 1 

Geophysics Division 

Pacifxc Missile Test Center 
Pt- HUgu, CA 93042 

31. Dr. Michael J- Krause 1 

AFGL/LIS 



Hanscom AFB, MA 01713 



55 



1 



32. HAJ Bob Bright 
AWS /DOOE 

Scott AFB, IL 62225 



33. HAJ Ed Kolczynski 1 

AWS/SYX 

Scott AFB, IL 62225 



34. Joel S. Davis 1 

Defense Sciences Div. 

Science Applications, Inc. 

1010 Woodman Dr.. Suite 200 
Dayton, OH 45432 



35. Mr. L. Biberman 1 

Institute for Defense Analysis 

400 Army Navy Dr. 

Arlington, VA 22202 

36. Dr. Richard Gomez 1 

DELAS-EO-MO 

Atmospheric Sciences Laboratory 
White Sands, NM 88002 

37. Dr. R. Fenn 1 

Air Force Geophysics Laboratory 

Hanscomb AFB, HA 02173 

38. Mr. Glen Spaulding 1 

MAT 72 

Naval Material Command 
Washington, D.C. 20362 

39. Dr. Lowell Wilkens 1 

Naval Weapons Center 

China Lake, CA 93553 

40. Dr. Ed Monahan 1 

Department of Oceanography 

University College 
Gallway, IRELAND 

41. Mr. Walter Martin 1 

Code 470, Rm 307 

Office of Naval Research 
800 N Quincy St. 

Arlington, vA 22217 



42. Dr. Gloria Patton 1 

Office of Naval Research 
1030 E Green St. 

Pasadena, CA 91106 



43. Dr. Warren Denner 1 

Science Applications Inc. 

2999 Monterey-Salinas Hwy 
Monterey, CA 93940 



56 



44 



1 



. Br. George Hanssen 

Science Spplicationc Inc. 

PO Box 1303 

1710 Goodrich Dr. 

BcLean, VA 22102 



45. Or. Lou Goodman 1 

Code 481 

Office of Naval Research 
Physical Oceanography 
NSTL Station, MS 39529 



46. CDR S.G. Colgan 1 

Code 420 B 

Office of Naval Research 
800 N Quincy St. 

Arlington, VA 22217 



47. Dr. John A. Cooney 1 

Department of Physics and Atmospheric Science 
Drexel Oniversity 
Philadelphia, PA 19104 



48. Officer in Charge 1 

Attn : LT J.P. Gleason, OSN 
NOCD, NAP Sigonella 
Catania, SICILY 
FPO New York 09523 



57 




.-CTA.: 



Thesis 
G4575 
c. 1 



. > ^ T 

Gleason 

Single-station 
assessments of the 
synoptic-scale forcing 
on the marine atmos- 
pheric boundary layer. 



